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For a quantum many-body problem, effective Hamiltonians that give exact eigenvalues in reduced 
model space usually have different expressions, diagrams and evaluation rules from effective transi- 
tion operators that give exact transition matrix elements between effective eigenvectors in reduced 
model space. By modifying these diagrams slightly and considering the linked diagrams for all the 
terms of the same order, we find that the evaluation rules can be made the same for both effective 
Hamiltonian and effective transition operator diagrams, and in many cases it is possible to combine 
many diagrams into one modified diagram. We give the rules to evaluate these modified diagrams 
and show their validity. 



I. INTRODUCTION 



Effective Hamiltonian iJ ff and transition operators 
O e ff are commonly used in many ab initio many-body 
calculations such as nuclear, atomic and molecular 
systems, and in phenomenolo gica l calculations 0, 0, 
SSSfiSIIillllHEi 13111 of doped transitional 
metal ions, lanthanide and actinide ions. H c g is defined 
to be an operator acting on a restricted model space of 
handleable dimensions to give upon diagonalization the 
exact eigenvalues and model space eigenvectors. For a 
time-independent transition operator O, O c g may be in- 
troduced that gives the same matrix elements between 
the model space eigenvectors of H c s as the original op- 
erator O between the corresponding true eigenvectors of 
H. Calculations based on H e g and O e ff have many ad- 
vantages over variational and other direct calculations 
based on H and O, such as smaller bases, less calcula- 
tion effort, order by order approximations etc, and can 
also be used together with variational and other direct 
calculations [? ] A recently review can be found in 

H c g and O e ff are often constructed by time- 
independent many-body perturbation theory (MBPT) 
with order-by-order approximation, which can be rep- 
resented with Goldstone diagrams, in analogy to Feyn- 
man diagrams|l8j. Many expansions have been given, 
whose diagrammatic representations usually contain not 
only connected but also disconnected diagrams, which 
have the shortcoming of size inconsistency, much more 
computation efforts and tremendous number of diagrams. 
Nonetheless, There are common-used effective Hamilto- 
nians, one hermitian and the other nonhermitian are 
known to contain only connected diagrams. The rules 
to generate diagrams and to evaluate them are well- 
known. Factorization theorem is shown to be able to 
combine diagrams having the same set of vertexes and 
lines but different relative orderings of vertexes, and 
hence reduces the number of high-order diagrams. Q 
Compared to effective Hamiltonians, effective transition 
operators generally have different algebraic forms, much 
more complicated diagrammatic representations, much 



greater number of diagrams and different diagram eval- 
uation rules. The hermitian (canonical) effective tran- 
sition operator, which works together with the hermi- 
tian (canonical) effective Hamiltonian, has been pre- 
sented in detail by Hurtubise and co-workers in a se- 
ries of papers. QUE! Duan and Reid[l9l| constructed a 
simple connected nonhermitian O e g that works together 
with the connected nonhermitian effective Hamiltonian 
and showed how to construct a connected expansion. 
Since it is well-known that hermitian effective Hamil- 
tonian up to third order can be obtained from trivial 
symmetrization, 6] and this can also be shown for effec- 
tive transition operators (up to third order in V, the 
perturbation in Hamiltonian), and also for higher-order 
calculations usually coupled-cluster methods come into 
play, most researches required only limited order dia- 
gram calculations of nonhermitian effective Hamiltonian 
and nonhermitian effective transition operators. How- 
ever, there are two problems with the effective operator 
methods: too many diagrams, and the rules to calculate 
energy denominators for effective transition operators be- 
ing different from those for effective Hamiltonian. [l6| 



Our aim in this paper is to modify the diagram rep- 
resentation of many-body perturbation expansion of ef- 
fective Hamiltonian and effective transition operators, so 
that the energy denominator rules for effective transition 
operators and effective Hamiltonian are the same, and 
the number of diagrams is greatly reduced and becomes 
handleable for third order terms of nonhermitian effec- 
tive transition operators. In section II we present our di- 
agrammatic representation of many-body operators and 
the corresponding evaluation rules for perturbation dia- 
grams; in section III we illustrate the generalized factor- 
ization theorem for the effective transition operators; and 
in Section IV we show how to denote several diagrams of 
the same set of vertexes and lines with only one diagram 
to group and reduce the number of diagrams. 
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II. MODIFIED DIAGRAMS AND THE 
EVALUATION RULES 



We use the following algebraic representation for a gen- 
eral fermion n-body interaction T n (l, 2, ■ • ■ , n): 



Tn = n\2 t QiQ2---a„,/3i/3 2 ---/3„ a i 1 ' ' - a i„ a /3„ ' ' 'aft, (1) 

ai,~ ,a„,0i,--,0 n 

where the coefficients t aia2 ... an ^ 1 p 2 ...fj n are anti-symmetry under exchange of a pair of bra indexes or ket indexes. It 
can be seen for one and two body interactions, the coefficients are given as follows 

t a , p = (^1)1^(1)1/3(1)), (2) 
t aia2 , 01 p 2 = (a 1 (l)a 2 (2)|T 1 (l,2)|/3 1 (l)/3 2 (2))-(a 1 (l)a 2 (2)|T 1 (l,2)|/3 2 (l)/3 1 (2)), (3) 



where the indexes in bracket "()" means the indexes of 
the variables of the function. It is straightforward to 
give the coefficients for n-body (n > 3)interactions. The 
diagrammatic representation of T n is: 



A 



a\\a 2 



Tn 



(4) 



The diagram is invariant under n\ interchange of the n 
out-going or in-going lines. The number of symmetry op- 
erations is (n!) 2 , which contributes to the factor l/(n!) 2 
as appeared in Eq. ^ For general close-shell reference 
"vacuum" state where there are both particle and hole 
lines, the number of inequivalent diagrams for T n are 
(n + l) 2 , corresponding to the n + 1 possibilities for the 
in-going lines (0, 1, ■ ■ •, n holes) and out-going lines. Par- 
ticle line and hole line are inequivalent, and hence the 
symmetry factor changes accordingly for each diagram. 

We need to modify the diagrams slightly so that the 
evaluation of energy denominators for effective operator 
is the same as that for non-hermitian Hamiltonian and 
is simple in general. This simplification helps us using 
generalized factorization theorem and simplifies diagrams 
developed in the next two sections. The modifications are 
shown in Fig and Fig [2 

The general rules for evaluating the diagrams of 
non hermitian Hamiltonian have been summarized by 
BrandowQ and Lindgren and Morrison [l8j. The rules 
of evaluating the diagrams for second order perturbative 
expansion of effective operator have been given by Hur- 
tubise and Freed [2£j. With the new notation QJ-@, 
the revised more general diagram evaluation rules can 
be summarized as follows: 1) draw only one diagram [4j 





(al) 



(a2) 



(1)) 



FIG. 1: A two-body third order diagram of non hermitian 
Hamiltonian, where (a2) is the original diagrams which con- 
tributes to the effective interaction showed by diagram (al), 
(b) is our modified representation for this diagram, where 
out-/in- going lines are changed into dashed line and end at 
dashed interaction vertex V e ff- 



for any one T n operator; 2) express the evaluation re- 
sult in the form Eq. (JJJ and the evaluation becomes 
calculation of iai,a 2 ,---;/3i-ftv coefficients; 3) in coeffi- 
cient <ai,- -,a m ;/3i,- -,/3 m i includes a factor (m!) 2 /(ii!i 2 ! • ■ ■), 
where m ingoing lines and m outgoing lines are partic- 
ipated into sets with i\, i 2 , • ■ •, equivalent lines, with i 
lines being equivalent if and only if they all start at the 
same interaction (or no beginning interaction), and point 
to the same ending interaction (or no ending interaction) ; 
4) include an over-all sign factor of (-l) l+h M 

where I is 

the number of closed loops and h the number of downgo- 
ing or "hole" line segments; 5) For each vertex, multiple 
the antisymetrized matrix elements as given in Eq. © 
and ©. 6) for a n-vertex diagram (suppose the vertex 
are labeled increasing from the bottom as 1, 2, • ••, n), 
include the products of energy denominators calculated 
just above the vertex 1, 2, • • •, n, where energy denom- 
inators are calculated by energy sum of all down going 
line taking away the energy sum of all upgoing lines, as 
show in Fig. ^ 7) sum each upgoing line independently 
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over all particle states, and each downgoing line indepen- 
dently over all hole states. The exclusion- violating terms 
which arise from these independent summations must be 



J 



included. 3 

The energy denominator can be calculated for diagram 
Fig.l (a2) as 



(£ 2 (l)-£ 2 (2))(£i(l)-£i(2)) 
= [(e p + e q ) - (-e„ - e b + e t + e n + e p + e q )] [{e p + e q ) - (-e b + e t + e s + e n )] (5) 
= (e Q + (b - e* - en)(e b + e p + e q - e t - e„ - e s ), (6) 

where Ei(j) is the energy for diagram Fig. ^ai) calculated at the indicated line by the energies of up- going lines 
taking away the energies of down-going lines. Diagram Fig. ^ (a2) can be written as 

(Vc)tn,ba(Vb)sa,pq(V a )mb,st \ 



Fig.|Ha2) = ^ E <n<a q a p f-(2!) 2 /2!^ — 

mn P q \ abts a 



+ e b - e t ~ e n )(e b + e p + . 



(7) 



r 



where the minus sign comes from the numbers of core 
solid lines, 2, and close cycle, 1, adding up to odd number 
3. the factor (2!) 2 /2! is due to the fact that the effective 
operator is two-body with equivalent lines p and q. 

The calculation of energy denominator can be put in 
an more systematic way as shown in Fig. \T]p: the energy 
denominator for i th vertex, is just the negative of net 



outflow energy (E noo ) of a loop enclosing vertex 1, 2, • • •, 
i, i.e., 

-E^eiVM ■■■V i ) = -(E aoe (V 1 )+E noe (V 2 )+- ■ ■+E Bm (V i )). 

(8) 

Therefore, the energy denominator factor for the n— ver- 
tex diagram can be written as 



D 



(-1) 



(-1) 



^E^MiE^Vx) + E noc {V 2 )] ■ ■ ■ [E^iYi) + E noc (V 2 ) + ■■■ E aoe (V n ^)\. 



r 



(9) 
(10) 




(al) 



FIG. 2: A second order (in V) diagram of effective transition 
operator O c ff, where V a and Vb are interactions contained in 
original Hamiltonian and O is the original transition operator. 



For Hermitian effective transition operators and 
Hamiltonians, the rules to evaluate diagrams have been 
given by Hurturbise and Freed 20J , which is different from 
Lindgren and Morrison's in energy denominators and 
contains algebraic factors. When we draw vertex of O c g 
at the same level as the vertex of O in perturbation di- 
agrams, as shown in Fig. Ufa), the rules of calculating 
denominators become exactly the same as those for the 
non-hermitian Hamiltonian. We modify Fig.^a) accord- 
ingly to give Fig. |2Jb) . In Fig. [2Jb) we save the free lines 
for the sake of calculating matrix elements. Note that 
the free lines should not be included in the calculation 
of denominators and the dashed lines should not be in- 
cluded in the calculation of matrix elements. Diagram 
Fig. |2fa) can be calculated as 



rmipq 



abst 



(V3) m b,st(Vl)nt,abO a 



s, P q 



+ £b - £t 



(11) 



It can been seen that the rules to calculate the denomi- are given in Eq. 
nators are the same for the two types of diagrams, which 
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FIG. 3: An example connected diagram with disconnected 
part at the bottom and on the top to illustrate factorisation 
theorem. 



III. GENERALIZED FACTORIZATION 
THEOREM FOR THE MODIFIED DIAGRAMS 

When several diagrams have the same set of vertexes 
and the same line directions and line types but different 
relative orderings of vertexes, they have the same matrix 
elements, factors and signs but different denominators. 
If all possible relative orderings are included, the results 
can be expressed by means of a single diagram with the 
denominators determined independently for each part of 
the diagrams 0. This is the factorization theorem. Lind- 
gren gave a illustration and proof H3- However it was 
not trivial to show that this theorem also holds for Hur- 
tubise and Freed's rules of evaluation of denominators for 
effective interaction diagrams. Here, the rules to evaluate 
denominators for our modified diagrams are the same for 



both effective Hamiltonian and effective operator, show- 
ing trivially that factorization theorem also holds for ef- 
fective operator. 

A diagram to illustrate this theorem is given in Fig. ® 
where V e and Vf are two disconnected parts at the bot- 
tom, and V a and Vb V c are two disconnected parts on the 
top. If one lowers V e relative to Vf, the vertexes and all 
the lines (direction and types) do not change, but the de- 
nominator from the lowest two interactions change from 
E noe {V f )[E noe (V e ) + E noc (V f )} to E noe (V e )[E noe (V f ) + 
E n oe(V e )]- By adding those two diagrams together, 
one gets a denominator E noo (Vf) E aoe (V e ), which is 
the same as calculating the lowest two disconnected 
interaction separately. The denominator for Vb is 
-E noc (VaV c V d V e V f ) = E noe (V b ), which follows from the 
fact that all -E^oc's for all vertexes add up to zero. Sim- 
ilarly, the denominator from the highest three interac- 
tions is E noc (V b )E noc (V b V a )E noc (V h V a V c ) for the order- 
ing showed in Fig. |5J There are two other diagrams 
with different orderings of V a relative to VbV c . By adding 
all the three diagrams up, one gets, similar to the case 
of the lowest two interactions, a denominator from the 
highest three interactions E noe {V a )E noe {Vb)[E noe {VbV c )]. 
The total denominator calculated this way for the sum 
of those diagrams with all the possible orderings is 
E noc (V a )E noc (V b )E noc {VbV c )E noc (V e )E noc {V f ). If one 
uses factorization theorem for both the upper discon- 
nected part and the lower disconnected part, one gets ex- 
actly the same denominator much more straightforward. 
For completeness, we give the result for this diagram as 
follows: 



Fig.® 



rnc,sn (Vb)e,u(Vc)ud,te{Vd)ab,cd(Ve)s,a(Vf)tb 

a m an\{ ) 2^ E noc (V a )E noc (Vb)E noc (V b V c )E noc (V e )E noc (V f ) 

abcde,stu 

E+ I ( (Va)rnc,ns(Vb)e.u(Vc)ud,et(Vd)ab.cd(Ve)s.a(Vf)tb 
a m a n \[ ) 2^ E noc (V a )E noc (Vb)E noc (V b V c )E noc (V e )E noc (V f ) 

abcde,stu 



\abcde,stu 



(V a ) (Vb)e.u(V c )ud,et(V d )ab,cd(Ve)s,a(Vf)tb 

(e c + e m - e„- e s )(e e - e„)(e„ + e d - e t - e m )(e e - e s ){e b - e t ) 
I 



(12) 
(13) 
(14) 



The sign (— l) 6 comes from the fact that there are five 
backwards solid internal lines and one loop (edbtue). If 
one switches the entering points of s and n at vertex V a , 
and e and t at vertex V Cl then one gets the equivalent 
diagram with three loops. The result for this diagram is 
given in 1)13(1. 



IV. DIAGRAMS WITH THE SAME SET OF 
VERTEXES AND LINE DIRECTIONS AND 
NOTATIONS 

The diagrams in Fig. 0] (b), (c), (d) and (e) have 
the same set of vertexes, line connections and directions. 
One can simply denote the four diagrams as diagram (a), 
where i and j can be both core lines and virtual lines. 
Hence the matrix elements, the symmetry factors and 
the E n0 cS have the same form. The denominators are 
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in,, f 







FIG. 4: An example to demonstrate the combination of sev- 
eral diagrams with the same set of vertexes, line connections 
and directions into one diagram. Here (b) has two particle- 
state internal lines, (c) and (d) each has one hole-state internal 
line and one particle-state internal line, (e) has two hole-state 
internal lines, and (a) is a diagram used to denote all the four 
diagrams (b) to (e). 



just different combinations of E n0 cS. One cannot com- 
bine the denominators together the same way as in the 
factorisation theorem since the orbital types for different 
diagrams are different. However, this does not prevent 
us from using (a) to represent all the four diagrams as 
long as one is cautious about the relations between the 
relative orderings of vertexes and the internal line types 
in evaluating diagram (a). 



Fig- El = J2a+a n ^r?i(i,j) 

run \ i,j 



mih])Eno C {V)E noc {W) 



where rji(i,j) = 1,-1,-1,1 and r]2{i,j) = — 1,1,1, — 1 
are the sign contributed by core orbitals and denomina- 
tors respectively. Note that factorisation theorem has 
been applied in the evaluation of diagrams (c) and (d). 



Denote the 
quentially V 2 
V n as (V n V n 
tial V 1 ,V 2 ,---,V n ,V 1 
(V 1 V n ---V 2 V 1 ), and 
then B as [AB]. 
Fig. 0t is uniquely 
gram Fig. 0J>e are 
DWV],(VDW)[VWD 



(e„ - ej)(e m - e t ) j ' 



(15) 



line entering V\ and then passing se- 
•• 5 ^7j_i and finally going out from 
•l • • • V^V^.), the loop passing sequen- 
and pointing from V± to V 2 as 
the ordering that A is higher 
It can be seen that diagram 
specified by (VDW) and dia- 
(VDW)[VDW],(VDW)[DVW + 
WVD] and (VDW)[WDV]. As 



has been argued above, one only needs to draw one dia- 
gram to work out the values of all those similar diagrams 
by using the denominator rule Eg. 1101 and other standard 
rules for matrix elements and factors. The application of 
these new notations to one and two-photon transitions 
can greatly reduce the number of diagrams. This will be 
presented elsewhere latter. 

V. CONCLUSION 



We have modified the diagrams for non-hermitian ef- 
fective Hamiltonians and effective operators so that the 
rules to calculate the denominators are the same for the 
two types of diagrams. Expressed with these modified di- 
agrams, proof of factorization theorem for effective tran- 
sition operators becomes trivial. By denoting the class 
of diagrams of the same set of vertexes, line connections 
and directions with one diagram, the number of diagrams 
can be greatly reduced and hence making the high order 
contributions to effective operators more handleable. 
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